Shear-melting of a hexagonal columnar crystal by proliferation of dislocations 
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A hexagonal columnar crystal undergoes a shear-melting transition above a critical shear rate 
or stress. We combine the analysis of the shear-thinning regime below the melting with that of 
synchrotron X-ray scattering data under shear and propose the melting to be due to a proliferation 
of dislocations, whose density is determined by both techniques to vary as a power law of the shear 
rate with a 2/3 exponent, as expected for a creep model of crystalline solids. Moreover, our data 
suggest the existence under shear of a line hexatic phase, between the columnar crystal and the 
liquid phase. 
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Understanding the behavior of complex fluids under 
flow is essential for their processing and technical use . 
The relevance of non-linear rheology for a wide range of 
applications has motivated many fundamental studies, 
both theoretical and experimental. However, the rheol- 
ogy of structured fluids is not well understood yet, despite 
robust experimental facts such as shear-banding associ- 
ated with marked modifications of the structure or tex- 
ture. This very general behavior is observed in a large 
variety of materials including amorphous, crystalline or 
liquid crystalline systems. Among liquid crystals, two- 
dimensional (2D) solids, or hexagonal columnar crystals, 
have been the subject of very few studies, but these ma- 
terials should deserve a particular attention because of 
both applicative and fundamental grounds. Their close 
analogy, in terms of elasticity and phase diagrams, to 2D 
magnetic flux line lattices (FLL) in type II high-T c super- 
conductors is appealing. On the other hand, as colum- 
nar liquid crystals are extensively used for the synthesis 
of mesoporous materials 2] , understanding and in turn 
controlling their structure under flow may lead to some 
technological advancements. 

We use synchrotron small-angle X-ray scattering 
(SAXS) and rheology to investigate the behavior under 
shear of a soft hexagonal columnar crystal, which con- 
sists of oil tubes arranged on a triangular lattice in water. 
Rheological experiments show a shear-induced transition 
between two states of markedly different viscosities. The 
low viscosity high shear rate structure is a 2D liquid of 
tubes, thus indicating that, above a critical shear rate or 
shear stress, a shear can induce the melting of the long- 
range 2D order of the tubes Q- We exploit the analogy 
with FLL to derive a scenario for the melting of the soft 
hexagonal columnar crystal and propose the melting to 
be due to a proliferation of dislocations. We combine the 
analysis of the shear-thinning regime below the melting 
transition and that of the SAXS patterns under shear to 
extract the shear rate dependence of the density of dis- 
locations, p. Both techniques show that p increases with 
the shear rate as a power law with the same exponent 
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FIG. 1: 2D diffraction patterns for a shear rate (left) 7 = 
3.2 sec -1 and (right) 7 = 1490 sec -1 in (top) tangential con- 
figuration with incident beam parallel to the velocity V, (bot- 
tom) radial configuration with incident beam parallel to the 
velocity gradient. Scale bar is 0.2 nm" 1 . 



2/3, as expected for a creep model of crystalline solids. 
The two independent determinations are fully consistent 
and hint at a dislocation-mediated melting of the hexag- 
onal columnar crystal under shear, in noticeable analogy 
with theoretical predictions for FLL £| . Our data more- 
over suggest the occurrence of an intermediate hexatic 
phase under shear. 

The experimental system is a lyotropic hexagonal 
phase consisting of infinitely long oil tubes immersed 
in water. The samples are composed of a quaternary 
mixture of sodium dodecylsulfate (SDS), pentanol, cy- 
clohexane and brine with a NaCl concentration of 0.4M. 
We use a composition in weight percent of 9.7% SDS, 
24.1% brine, 61.5% cyclohexane and 4.7% pentanol which 
yields, at rest, a structure consisting of oil tubes of radius 
15 nm, coated with a surfactant monolayer and arranged 
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FIG. 2: Flow curve, stress versus shear rate, in the low shear 
rate regime. Curves are labeled by temperature. 



on a triangular array with a lattice parameter clq = 33 nm 
Rheology experiments are performed in a Couette 
geometry with a stress-controlled Paar Physica UDS 200 
rhcomctcr. The flow curve, stress a vs shear rate 7, of 
a soft hexagonal columnar crystal, has been described in 
detail in Refs.0, 0. At low shear rates (7 < %) the 
system behaves as a power law shear-thinning material 
while at high shear rates, the stress varies affinely with 
the shear rate. These two regimes correspond to two 
stable branches of stationary states, for which data ob- 
tained by imposing either a or 7 exactly superpose. The 
transition from the lower branch to the higher branch oc- 
curs through a stable hysteretic loop in a stress-controlled 
experiment. To characterize the structure of the ma- 
terial under flow, SAXS experiments under controlled 
shear rate in Couette geometry are performed on the 
ID-2 beamline at the ESRF, Grenoble, France 1\. The 
structure of the sample in the two stable branches has 
been previously characterized Examples of the 2D 
patterns obtained in both tangential and radial configu- 
rations are given in Fig. ^ At low shear, the hexagonal 
phase exhibits a polycrystalline texture with the tubes 
progressively aligned along the velocity direction as 7 in- 
creases, although a state with the tubes perfectly aligned 
along the flow is never reached. The angular width of the 
arcs (Fig. bottom) is used as a measure of the align- 
ment and indicates a sample mosaic spread, A9, which 
decreases monotonically from 26deg to 15deg when the 
shear rate increases from 1 to 100 sec -1 . By contrast, at 
high shear rate, the long-range two-dimensional order of 
the tubes melts, leading to a 2D liquid of tubes strongly 
aligned along the flow (A8 ~ 8deg). Melting transitions 
under shear have been frequently observed for three di- 
mensional colloidal crystals and are generally believed 
to be due to the existence of a periodic potential in the 
direction of flow || . Similar arguments can certainly not 
be invoked for a 2D solid, for which a crystalline array of 
tubes perfectly aligned along the velocity should flow at 
low viscosity. A novel approach is thus needed. 



To better understand the physical mechanism for the 
shear-melting, the behavior of the material under mod- 
erate shear is analyzed within the framework of work 
hardening of crystalline solids. The parallel between the 
behavior under shear of liquid crystal phases and that of 
metals has been proved successful in the case of lamel- 
lar systems |l0j . In the formal theory of work hardening, 
the stress required to move a dislocation through a region 
having a density of dislocation p can be derived from gen- 
eral dimensional arguments 11]. With the assumption 
that the applied stress, a, is relaxed by the dislocations, 
a should depend solely on three parameters, the density 
of dislocation, p, the Burger vector of the dislocations, b 
and the shear modulus of the material, G. The pure num- 
ber a/G must therefore be a function of the pure number 
bp 1 / 2 . On the other hand, the force which a exerts on a 
unit length of a dislocation, the Peach and Koehler (PK) 
force, is ba. This force is balanced by a resisting force 
arising from the line tension, which is proportional to b 2 . 
In a stationary state, ba ~ b 2 . Thus, a varies linearly 
with b and a relation between the stress and the density 
of dislocation follows: 

a = KGbp 1/2 (1) 

where if is a numerical factor. Equation JQ) is well 
obeyed by crystalline solids, with the prefactor K of 
the order of 2 for metals |l2(. For hexagonal colum- 
nar phases, there are three generic types of dislocation 
[IsT ] : a screw dislocation, a (longitudinal) edge dislocation 
with tangent t parallel to the columns dislocation and a 
(transverse) edge dislocation with t perpendicular to the 
columns. The transverse edge dislocation is the most 
costly in energy because its creation requires the forma- 
tion of two end-caps and is thus presumably very rare 
in soft columnar crystals. By contrast, longitudinal edge 
and screw dislocations do not require the creation of end- 
caps and are certainly more easily created; these disloca- 
tions are expected to form loops with both longitudinal 
edge and screw components, as also assumed in Ref. jl4| 
for FLL. For these two types of dislocation 0], the line 
tension is proportional to b 2 . Thus, Eq. Q should hold 
and an increase of the density of dislocation with the ap- 
plied stress is predicted. Moreover, whenever dislocation 
motion is the dominant plastic deformation mechanism, 
one observes a constant shear rate regime usually de- 
scribed by Orowan's geometrical relation, 7 = pbv |l5l |. 
where v, the average velocity of the dislocations, is pro- 
portional to the applied stress a: v = Ma, with M a 
mobility. The shear rate dependence of both the stress 
and the density of dislocations can thus be extracted from 
the combination of Orowan's relation and Eq. One 
obtains: a = i^) 1 ' 3 x 7 1 / 3 and p = (bGM)~ 2 / 3 x 7 2 / 3 . 
We are able to measure independently the shear rate de- 
pendences of both a and p. As we shall see below, we 
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find p ~ 7 2 ™ and a ~ 7™ with n = 1/3, in remarkable 
agreement with the simple dimensional theory described 
above. 

Rheology measurements give directly the shear- 
thinning behavior of the sample. As can be seen in 
Fig. |21 in a large range of temperature, from 15 to 
30 °C, and over more than one order of magnitude for 
7, we find a ~ 7" with n = 0.32 ± 0.03, which implies 
p ~ ^,o.64±o.06_ rpj^ gjjggj. ra ^- e var j a tion of the density 

of dislocation can be confirmed by independent measure- 
ments. Thanks to the Synchrotron high resolution, quan- 
titative information on p can indeed be obtained from the 
width of the diffraction peaks. The peaks are not reso- 
lution limited and the full width at half height (FWHH) 
of the first order diffraction peak, Aq, increases with 7, 
while its position remains unchanged. Figure shows 
that the FWHH ranges from 10 _2 nm _1 at rest up to 
2.2 x lO^nm" 1 for 7 = 100 sec" 1 . We note that this 
curve is fully reversible and in particular the zero shear 
value of Aq is systematically recovered when the shear is 
stopped. We find that the variation of Aq can be very 
well fit by the sum of a constant and a power law of 
the shear rate: Aq = Aqo + 5(7)1 where the constant 
Aqo is related to the size of the crystallites at rest. We 
obtain Aq = 0.9 x lO^nm" 1 and 3(7) = B^™ 1 , with 
B = 0.3 x lO^nm" 1 and m = 0.33 ±0.06. Finite size ef- 
fects can account for the shear rate dependence broaden- 
ing of the diffraction peaks. The standard Laue-Scherrer 
relation relates indeed the width of the Bragg peaks to 
£, a translational correlation length: 3(7) = 27r/£. In the 
case of a crystalline solid, £ can be identified as the aver- 
age distance between dislocations and thus p = £" 2 . 
From the power law variation of g, one predicts the den- 
sity of dislocation p to scale as 7 2 " 1 = ^°- 66 ± 08 . Re- 
markably, the same power law variation of the density of 
dislocation with the shear rate with an exponent 2/3 is 
deduced from the SAXS data and from the rheological re- 
sults, thus providing convincing support to our scenario 
for plastic deformation. 

The variation of N = £/<Zo; the translational corre- 
lation length (normalized by the lattice parameter hq), 
with the shear rate, below and above the melting tran- 
sition, are reported in Fig. 0] A continuous decrease of 
N with 7, from 80 to 4, is measured. Our results clearly 
show two regimes characterized on a log-log plot by a 
change of slope at the shear rate j c . The slope is equal to 
— 1/3 below 7 C and crosses over a more abrupt variation, 
where a slope of —0.64 ± 0.06 is measured. Although the 
cross-over at the melting transition is smooth, it is asso- 
ciated with a dramatic drop of viscosity and signs unam- 
biguously the melting transition. At 20 °C, temperature 
at which the SAXS experiments are performed, both the 
rheology and the SAXS experiments show that the shear 
melting occurs for j c ~ 150 sec" 1 . It is associated with a 
critical normalized correlation length N c = |^ ~ 14 (Fig. 
0}, which is interestingly of the same order of magnitude 
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FIG. 3: Shear rate dependence of the full width at half height 
of the first diffraction peak of the sample. The solid line is 
the best fit of the data (see text). Insert: First diffraction 
peak for a sample at rest (squares), submitted to a shear rate 
7 = 5.8 sec" 1 (circles) and 7 = 29 sec" 1 (triangles). 
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FIG. 4: Shear rate dependence of the normalized translational 
correlation length, as measured from the broadening of the 
first order Bragg peak. The two solid lines have slopes of 
-1/3 and -0.64. 

as the correlation length measured at the melting of a 
2D colloidal crystal 17]. Additionally, for b = ao, which 
are presumably the most frequent dislocations, the nu- 
merical prefactor K in Eq. is equal to <7 '^ c , where 
er c ~ 23 Pa (inset Fig. ^| is the critical stress at the melt- 
ing transition and G ~ 250 Pa is the shear modulus pre- 
viously measured 0]. One thus obtains K ~ 1.3, in 
excellent a gree ment with the value experimentally found 
for metals [12J. 

Both the shear thinning regime and the SAXS scat- 
tering data can be quantitatively and consistently inter- 
preted as resulting from an increase of the density of dis- 
locations with the shear rate. We therefore propose the 
shear melting be due to a proliferation of dislocations. 
Our out-of-equilibrium results are in remarkable anal- 
ogy with the shear melting of FLL Q , theoretically ad- 
dressed when the shear is perpendicular to the flux lines. 
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FIG. 5: (Right) 2D diffraction pattern of a sample submitted 
to a shear rate 7 = 390 sec -1 . The velocity is perpendicu- 
lar to the scattering plane. (Left) Azimuthal scan of the 2D 
diffraction pattern. The intensity I is averaged over an an- 
nulus of radius go = 0.215 nm -1 and width 0.05 nm" 1 and is 
corrected to take into account the curvature of the cell; once 
corrected, a six fold symmetry is recovered with all maxima 
of roughly equal intensity. 



In our case, however, the shear is oblique to the tubes; 
the non trivial coupling between the gradual alignment 
of the tubes and the increase of the density of dislocation 
with the shear remains an open issue. Our experimen- 
tal results compare also naturally with those for the well 
documented melting of 2D solids (lS|. However, because 
scattering experiments do not provide direct imaging, no 
information about the spatial distribution of the disloca- 
tions can be accessed, but only their average density can 
be extracted from the data. As opp osed to 2D solids, 
where direct imaging is possible 20], it is therefore 
delicate for hexagonal columnar crystals to unambigu- 
ously characterize the melting transition, and determine 
whether it is a one-stage first order transition or a two- 
stage transition with an intermediate line hexatic phase. 
Our data tend nevertheless to support the latter scenario. 
Indeed, fig. [S] shows the 2D pattern, obtained in tan- 
gential configuration for 7 = 390 sec -1 , that is above 
the melting transition. The pronounced six-fold angular 
modulation observed on this pattern is a signature of a 
line hexatic phase. Line hexatic phases have been pre- 
dicted [lj, |21j and only very recently observed by scat- 
tering techniques |22ll23| . Interestingly, our experimental 
results also suggest the existence under shear of an inter- 
mediate line hexatic phase, between the columnar crystal 
and the liquid phase. 

To conclude, we have shown that a mechanism of pro- 
liferation of dislocations can qualitatively account for the 
melting of a hexagonal columnar crystal under shear, 
which presumably occurs through a two-step process with 
an intermediate line hexatic phase. Our experimental re- 
sults compare well with the shear-melting of FLL and 
provide an original illustration of the deep analogy be- 
tween magnetic flux line lattices and columnar liquid 
crystals. Recent developments in the simulations of dis- 



location motions should render possible a numerical test 
of our main findings and should provide decisive infor- 
mation about the spatial organization of the dislocations 
in a steady-state regime. 
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